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We propose a hexagonal optical lattice system with spatial variations in the hopping matrix 
elements. Just like in the valley Hall effect in strained Graphene, for atoms near the Dirac points the 
variations in the hopping matrix elements can be described by a pseudo-magnetic field and result 
in the formation of Landau levels. We show that the pseudo-magnetic field leads to measurable 
experimental signatures in momentum resolved Bragg spectroscopy, Bloch oscillations, cyclotron 
motion, and quantization of in-situ densities. Our proposal can be realized by a slight modification 
of existing experiments. In contrast to previous methods, pseudo-magnetic fields are realized in a 
completely static system avoiding common heating effects and therefore opening the door to studying 
interaction effects in Landau levels with cold atoms. 

PACS numbers: 73.22.Pr, 73.43.-f, 67.85.-d, 37.10.Jk 


The Lorentz force, which acts on charged particles 
moving in a magnetic field, results in a number of fun¬ 
damental phenomena in condensed matter systems in¬ 
cluding the Hall effect in metals, Abrikosov lattices in 
superconductors, and the integer and fractional quan¬ 
tum Hall effects in ultra-pure two-dimensional electron 
gases. While phenomena, such as the quantized conduc¬ 
tance plateaus of the integer and fractional quantum Hall 
effects have been both observed experimentally and de¬ 
scribed theoretically [T], many properties, such as the 
non-abelian nature of excitations in the fractional quan¬ 
tum Hall effects Em, remain subjects of active research. 

These problems are difficult - they involve strongly in¬ 
teracting systems that are resistant to conventional theo¬ 
retical and numerical tools. The current theoretical state 
of the art involves comparisons of trial wave functions 
and numerical calculations on small systems using exact 
diagonalization and the density matrix renormalization 
group (DMRG) [3l|4]. Ultra cold atom experiments offer 
an alternative route, in which, potentially, the interplay 
of gauge fields, band structure, interactions, and disor¬ 
der can be studied by engineering and controlling these 
effects independently [5]. Moreover, by engineering these 
properties one could generate novel phases that have not 
yet been observed in condensed matter systems 

Various groups have recently experimentally demon¬ 
strated ‘synthetic gauge fields’ - methods for driving neu¬ 
tral atoms using laser beams in such a way that they be¬ 
have as if they were charged particles moving in a mag¬ 
netic field [EHI2- A number of effects, such as Abrikosov 
lattice formation [13], Hall deflection ED H H], and 
chiral currents [T6|, have been observed. However, an 
important limitation of these methods, that use either 
periodic lattice modulations or Raman transitions, seems 
to be significant heating of the atom clouds. For op¬ 
tical lattice experiments, this limits the timescale in 
which experiments can be performed to several tens of 


milliseconds ED, in contrast to experiments in static lat¬ 
tices which allow for several hundreds of milliseconds EH]. 
Further, experiments so far have been performed with 
noninteracting or weakly interacting atoms. An exten¬ 
sion of synthetic gauge fields experiments to the strongly 
interacting regime would require low heating rates in 
combination with low initial temperatures. 

In this Letter, we propose an alternative method for 
generating synthetic magnetic fields in ultracold atom 
systems that relies on a completely static optical lattice, 
and leads to the formation of relativistic Landau levels. 
Inspired by pseudo-magnetic fields observed in strained 
graphene EaEOj, molecular graphene [21] , and photonic 
systems [22], we propose a method for generating spa¬ 
tially varying hopping matrix elements in a hexagonal 
optical lattice. Starting from a standard configuration 
of three Gaussian laser beams intersecting at 120°, our 
method relies on simply displacing the beams [Fig. m- 
We show that pronounced Landau levels are generated 
close to the original Dirac points corresponding to al¬ 
most homogeneous magnetic fields in real space with 
opposite sign at the nonequivalent Dirac points. This 
leads to a host of observable phenomena, such as the val¬ 
ley Hall effect, quantization of in-situ densities, Landau- 
Zener effects in Bloch oscillations, and the emergence of 
gap structure in Bragg spectroscopy. Since our method 
relies on fully static optical lattices, it forms an attrac¬ 
tive starting point for investigating interaction effects in 
Landau levels [iiaEUlHI with ultracold atoms. 

Pseudo magnetic fields in optical lattices - Consider an 
optical lattice that is created by the intersection of three 
blue-detuned laser beams at 120° angles ESEOj. The 
resulting honeycomb potential has the form 
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FIG. 1. (a) Schematic of the setup: the optical lattice is produced by three Gaussian laser beams intersecting at 120° with 
offset d - lattice orientation is depicted in the upper right. The cross marks the center of the harmonic trap, and the star 
marks the position associated with the displaced Dirac cones depicted in (b). (b) Schematic of the displacement of the Dirac 
cones in momentum space associated with stretching type 1 bond. First Brillouin zone is indicated with the dashed line, (c) 
Pseudo-vector potential A as a function of position, (d) Pseudo-magnetic field 5 = V x ^ as a function of position. The 
hexagon in (c) and (d) marks the sample area with 80% uniformity in the pseudo-magnetic field. 


where km is the wave vector of the m-th laser beam and 
/^(r) is it’s position dependent intensity that accounts 
for the Gaussian nature of the beams. We consider the 
case that the lattice is sufficiently deep so that the tight 
binding model is applicable. If all three beams have 
identical intensities the resulting band structure (in the 
lowest band) is identical to that of graphene with two 
degeneracy points in the first Brillouin zone described 
by the Dirac equation. Changing the intensity of one of 
the beams corresponds to applying a uniform strain to 
graphene m, which can be captured by modifying the 
hopping matrix elements in the tight binding model and 
results in the shift of the Dirac cones in the Brillouin zone 
[see Fig. Ef]- This shift can be encoded in terms of the 
Dirac equation by adding a vector gauge field A: 

i^Dirac = fiVf {-id^ + A^) (J^ , (2) 

where Vf is the group velocity near the Dirac point, and 
are the Pauli matrices. Allowing for the intensities of 
all three beams to vary in space, we can obtain a non- 
uniform A{f). 

The tight binding model with non-uniform hopping 
matrix elements on the honeycomb lattice is [32] 

H = — tija\aj + V^a|a^, (3) 

(bi) * 

where tij are the hopping matrix elements, Vi = 
is the on-site potential that is a combination 
of the trap potential and the anti-trapping effect of the 
blue lattice uj‘1^ = cj^rap “ ^anti-trap with TTiat the atomic 
mass, and a\ and are the creation and annihilation 
operators. There are three types of hopping matrix ele¬ 
ments tij associated with the three hopping directions in 
a honeycomb lattice. We label these as tij = tu{r), where 
(g + g )/2 and u is the index of the laser beam that 
is perpendicular to the vector fi — fj (see Fig. [^). Us- 
ing this notation, and making the approximation that 


ti(r) t 2 {r) ts{f) and all are slowly varying functions 
of r, we find 

/ A^{f} \ ^ , V3 / 2t3{f) - ti{f) - tiif) \ 

\Ay{f)) 2\to{f)\ y3(i2(r-)-U(r)) T 

Here, A is the wavelength of the laser, to(r) = (G(^) + 
^ 2 (^) + and ± corresponds to the gauge field at 

the two distinct Dirac points. We note that this form of 
the vector gauge is identical to that obtained in strained 
graphene, but the origin of the variation in the hopping 
matrix elements is different - here it is induced by spa¬ 
tial variation of the lattice depth as opposed to strain in 
graphene [33l[34] , 

We connect the tight binding parameters tu{r) and 
V{f) to the laser light intensity lui'k) at point r using a 
simple analytical model that accurately captures a pre¬ 
cise numerical calculation [35|. Our model El assumes 
that the connection is completely local [valid when Iu{r) 
varies slowly in space on the length-scale of a unit cell] 
and that the laser beam intensities are close to isotropic 
[/i(r) ~ precisely the assumptions re¬ 

quired for the validity of the pseudo-magnetic field de¬ 
scription. 

A prescription for a uniform pseudo-magnetic field - 
To introduce spatial variations of beam intensities, we 
propose to use three Gaussian beams with the same in¬ 
tensity and beam waist and to align them such that the 
beam axis form an equilateral triangle Fig. [^. This 
choice results in a symmetric gauge field depicted in 
Fig. [5, and a nearly uniform pseudo-magnetic field de¬ 
picted in Fig. p. We note that this prescription does not 
introduce an offset between A and B sub-lattices. We can 
tune the strength and uniformity of the pseudo-magnetic 
field by varying the beam waist size wq and the displace¬ 
ment parameter d [Fig. [^]. In order to probe uniform 
phases of matter, such as the quantized Hall effects, one 
typically requires a uniform magnetic field over the sam¬ 
ple area. We have tabulated the optimal choice of wq and 
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FIG. 2. (a) Local Density of States as a function of energy and position in the trap for various trap frequencies, entrap = 
11.7(27rHz) - trap cancels the anti-trapping potential of lattice beams. Cc;trap = 16.4(27rHz) - bending of distinct Landau. 
Co’trap = 21.1 (27rHz) - strong smearing of Landau levels, (b) Density of a fermionic fluid as a function of position in the trap 
showing incompressible plateaus [corresponding chemical potentials are indicated with dashed lines in (a)]. 


d for a range of sample sizes Rq that ensure a uniform 
pseudo-magnetic field m- 

The maximum pseudo-magnetic field is limited be¬ 
cause the description in terms of the Dirac equation 
with A breaks down when the displacement of the Dirac 
cones becomes comparable to the the size of the Brillouin 
zone \A\ ^ 1/A. For the symmetric gauge choice \A\ 
varies linearly across the sample, and hence the max¬ 
imal pseudo-magnetic field is inversely proportional to 
the sample size. From Table II of the supplement we find 
^max = V X ^4 ~ 2.7/XRq. 

Proposed experimental setup - To make a concrete con¬ 
nection to experiment, we focus on a particular realiza¬ 
tion similar to the one in Ref. [30], which we will use 
throughout the remainder of the letter. We consider ^^Rb 
atoms in an optical lattice with A = 700nm, and hence 
a recoil energy of Eji = 2rAT?^ h x 4.685kHz. 
We choose the laser beam intensity such that the po¬ 
tential is 4:Er (per beam, in the beam center) which en¬ 
sures that the lattice is sufficiently deep to make the tight 
binding model applicable but sufficiently shallow that the 
hopping timescales are acceptably fast [Uj/h ~ 868Hz]. 
Finally, we focus on the case Rq = 23/im, which cor¬ 
responds to a hexagonal sample with 15,000 sites. For 
this setup, the maximum uniform pseudo-magnetic field 
of B = 0.17/im“^ [Fig. [^] is obtained by setting wq = 
150/im and d = Abiam. 

Landau levels in a harmonie trap - The effective trap¬ 
ping frequency cjeff must be large enough to ensure that 
the atoms are confined, but not so large as to smear out 
the discrete Landau levels. Without a trapping poten¬ 
tial, we expect to see ‘relativistic’ Landau levels at en¬ 
ergies En = sign(n) hvF^J‘2\n\B^ where vp = Xtij/{\/3h) 
and n is an integer [Ml EH- The corresponding Lan¬ 
dau level wave functions have a length scale set by the 
cyclotron radius Tc = {2n + f)/B. The length scale 

for the ground state of the trapped lattice system is 
^SHO = [2tijA^/(9mRb^^eff)] hence the ground state 
wave function will be confined if Asho ^ which sets 
a lower bound on entrap- On the other hand, a Landau 


level will avoid being smeared out if the shift of the trap 
potential on the scale of Tc is small compared to the spac¬ 
ing to the next Landau level ^ ^|n|+i ~R\n\^ 

which sets an upper bound on cjeff- Putting these consid¬ 
erations together, we find a window for observing Landau 
Levels 

‘2X^tij ^ 2 < 7.2tij f ^/\n + 1 | — 

9mRbRo rriRbRo^^"^^^^ \ ^ ^ / 

where in the last expression we have used the relation be¬ 
tween Rq and B. For our realization 0.2Hz < cjeff/^Tr < 
{91Hz, 12Hz} for n = 0, ±1. 

To visualize the Landau Levels in a trapped system 
we plot r), the local density of states (LDOS), as a 
function of position and energy [Fig. [^]. In the absence 
of a trap, we expect to see sharp peaks in the LDOS as a 
function of energy that are independent of position. For 
|VR(r)| \En — En-i\/rc^ the LDOS is simply shifted 
p{E,r) ^ p{E — R(r),0) while for large |VR(r)|, this 
picture breaks down and the peaks in p{E,r) become 
smeared. 

Experimental signatures - One of the defining char¬ 
acteristics of the quantum Hall effect is the formation 
of incompressible plateaus associated with the Landau 
levels. We propose that these plateaus can be observed 
with ultra cold fermionic atoms by measuring the atom 
density as a function of distance away from the center of 
the trap [Fig. |^]. In principle, it should be possible to 
observe multiple plateaus in a “wedding cake” structure. 
Due to the limitations on experimentally accessible opti¬ 
cal lattice sizes we expect to see only one or two plateaus. 

Another important limitation of the present generation 
of experiments with fermions in optical lattices is that the 
lowest temperatures that one can achieve T ^ f).btij [38l 
l39] are similar to the spacing of the Landau levels that 
we predict for the smallest samples m- An attractive 
alternative is to use bosonic atoms that can be Bose- 
Einstein condensed (BEG) in the lattice. 

To detect the pseudo-magnetic field we propose ob¬ 
serving cyclotron motion of a cloud of condensed atoms. 
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FIG. 3. (a) Snap shots of the in-situ density of an atom cloud undergoing cyclotron motion, where A^o is the total number of 
atoms in the cloud. Trajectory of the center of mass of the cloud is superimposed on top (red dots - tilt ajh — 0.84kHz//im, 
blue - no tilt, green - current frame), (b)-(c) Momentum resolved Bragg-spectroscopy: transition rate as a function of frequency 
V — ujjl'K and momenta kx and ky. (b) Slice at hxed ky=0. (c) Slices at hxed u (as indicated by the dashed white lines in (b)). 
The location of the Dirac cone is indicated with dashed magenta line in (b) and (c). 


Starting with the cloud in the ground state in the cen¬ 
ter of the trap, we remove the trap (set cjeff = 0) and 
accelerate the cloud towards one of the Dirac points by 
applying a tilt g • r (g is a unit vector indi¬ 

cating tilt direction, and a is the tilt ‘force’). When the 
momenta of the atoms approaches the Dirac point, we 
remove the tilt and observe the evolution of the cloud 
in situ. As the atom cloud is moving with a momentum 
close to the Dirac momentum, it is well described by the 
Dirac equation, and hence its trajectory is curved by the 
pseudo-magnetic field [Fig. [^]. The effective Lorentz 
force changes sign if the cloud is accelerated towards a 
non-equivalent Dirac point. 

In order to directly observe the Landau levels, we pro¬ 
pose to use Bragg spectroscopy. In this setup, a BEG that 
is originally prepared in the ground state 'ipo is trans¬ 
ferred to excited states -0^ using a Bragg spectroscopy 
setup [m ED] described by the perturbation potential 
V{r,t) = Vi cos(^ • f)cos{ujt). In Fig. |^b), (c) we plot 
the transition rate for a range of k and uj. As a func¬ 
tion of cj, we see clear maxima that correspond to the 
various Landau levels. Moreover, we can extract the mo¬ 
mentum structure of the Landau level wave functions by 
varying k (i.e. performing momentum resolved Bragg 
spectroscopy). Without the complexity of a Bragg spec¬ 
troscopy setup, Bloch-Zener spectroscopy mi offers an 
alternative that can detect the separation between the 
n = 0 and n = ±1 Landau levels HZ]. 

Observing interaction effects- Pseudo-Landau levels 
are particularly susceptible to the effects of interactions 
as the kinetic energy of the atoms within a single pseudo- 
Landau level is quenched. Previous theoretical inves¬ 
tigations showed that the interplay of short and long 


range interactions in strained hexagonal lattices can drive 
the formation of a number of exotic phases IS H IH- 
1^ . Specifically, at 1/2-filling of the zeroth Landau 
level, long range interactions prefer the formation of 
the quantum anomalous Hall effect (that spontaneously 
breaks time reversal symmetry) [2S Ell EM2H1 while 
short range interactions prefer the formation of charge- 
density-wave phases [2S EH] • Detailed numerical analy¬ 
sis of the 1/2 filling case with on-site interactions alone 
suggests the formation of an exotic mixture of ferro- 
and anti-ferromagnetism [7]. At 2/3 filling of the ze¬ 
roth Landau level, tuning the ratio of long and short 
range interactions is predicted to drive the formation of 
a triplet superfluid, a fractional topological insulator and 
the 2/3 fractional quantum Hall effect [6]. In compari¬ 
son to natural graphene, our proposal gives the advan¬ 
tages of (1) better control over the strain patterns and 
(2) control over both short [5] and long-rang interactions 
(using methods based on Rydberg atoms [42], dipolar 
atoms [43], or dipolar molecules [44]) and (3) control over 
filling factors and (4) the availability of low disorder po¬ 
tentials, making it particularly promising for realizing 
these exotic phases. 

Outlook- We have proposed a scheme for generating a 
pseudo-magnetic field in ultra cold atom systems that 
relies on spatial variations of the hopping matrix in 
analogy to the case of strained graphene. The typical 
timescales for conventional synthetic gauge field experi¬ 
ments in optical lattices is several tens of milliseconds m 
mostly limited by drive induced heating. Our approach 
could extend this to several hundreds of milliseconds as 
typical in static lattices [18], thus improving timescales 
for investigating the interplay of pseudo-magnetic fields 
and interactions. Generally, our work could establish a 
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new, very explicit link between solid state physics and 
ultracold quantum gases connecting ongoing experimen¬ 
tal work in strained graphene and novel, strained optical 
lattices. 
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